Helstrom's Theory on Quantum Binary Decision 

Revisited 



Gianfranco Cariolaro and Alberto Vigato 
Department of Information Engineering (DEI), University of Padova, via G. Gradenigo 6/B - Padova, PD 1-35131, Italy 

Email: {cariolar,alvigato} @dei. unipd.it 



o 

> 

O 

(N 



C/5 



> 
00 
00 

m 
in 

d 



X 



Abstract — For a binary system specified by the density opera- 
tors po and pi and by the prior probabilities qo and 51, Helstrom's 
theory permits the evaluation of the optimal measurement 
operators and of the corresponding maximum correct detection 
probability. The theory is based on the eigendecomposition (EID) 
of the operator, given by the difference of the weighted density 
operators, namely D — q^pi — qopo- In general, this EID is 
obtained explicitly only with pure states, whereas with mixed 
states it must be carried out numerically. In this letter we show 
that the same evaluation can be performed on the basis of a 
modified version of the Gram matrix. The advantage is due 
to the fact that the outer products of density operators are 
replaced by inner product, with a considerable dimensionality 
reduction. At the limit, in quantum optical communications the 
density operators have infinite dimensions, whereas the inner 
products are simply scalar quantities. The Gram matrix approach 
permits the explicit (not numerical) evaluation of a binary system 
performance in cases not previously developed. 

I. Introduction 

In 1970 Helstrom et. al. Q] evaluated the optimal per- 
formance (minimum error probability) of a binary quantum 
system, both with pure and mixed states. Till now, this result 
represents the only explicit solution of quantum optimization. 
In fact, for multilevel quantum systems the optimal solution 
is known in very few particular cases. 

The formulation of Helstrom's theory is the following. Let 
be po an d p\ the density operators of the given binary quantum 
system and let qo and qi be the corresponding prior probabili- 
ties. Then, form the decision operator as D — q\p\ — qopo and 
consider its eigendecompositon (EID), which has the general 
form 

D = QtPt ~ goPo = y^VkMiVk] (1) 
k 

where are the eigenvalues and \rjk) the corresponding 
eigenvectors. Finally, the correct detection probability is given 
by 

P c = qo+ ^2 (2) 

»?fc>0 

and the optimal measurement operators by 

nx=^kX%|, n =X>fcX%l- (3) 

Vk>0 i) k <Q 

However, the explicit evaluation poses severe difficulties when 
the dimensions of the underlying Hilbert space is larger than 
two 0. 



Here, we reconsider Helstrom's theory presenting a method 
based on a "modified" version of the Gram matrix, which 
we call skew Gram matrix. The net difference is that, while 
the operator D is expressed in terms of outer products (see 
dlj) and has the dimension of the Hilbert space, the Gram 
matrix is defined in terms of inner products. The advantage 
of the replacement of outer products by inner product is a 
considerable dimensionality reduction. At the limit, in quan- 
tum optical communications the density operators have infinite 
dimensions, as stated by Glauber's theory on the representation 
of quantum laser radiation field [3], whereas the inner products 
are simply scalar quantities. 

This new approach allows to find explicit results for the 
optimal detection, not only for the case of pure states (rank 1 
density operators), as found in 0J, &I, but also with mixed 
states. The explicit results are particularly simple when the 
quantum states have the geometrical uniform symmetry (GUS). 

II. Binary detection based on the skew Gram 

MATRIX 

Let ki = rank(pi) and assume that ko + ki < n, where n is 
the dimension of the underlying Hilbert space. Let 7, be the 
nx ki factors of p' i = qipu so that 7o7o = Pq and 7171 = p[ 
(* denotes the conjugate transpose). Then, the state matrix is 
given by r = [70, 71] and the Gram matrix by 

7o7o 7o7i 1 [Goo G01 
7i7o 7*7i J |_ G10 Gn 
where the blocks G^ = 7 t *7j has dimension ki x kj. The skew 
Gram matrix (SGM) is obtained by changing the sign of the 
first block row, namely 

-Goo — G01 
G10 G11 

Note that the ordinary Gram matrix G is Hermitian and 
positive semidefinite (PSD), and hence Goo, G11 are Hermitian 
and PSD, while G w = Gq V The SGM is neither Hermitian 
nor PSD. Note also that, given the density operator p' i7 
the factors 7,; are not unique, but the consequence of their 
multiplicity is irrelevant for quantum detection 0. 

Let us first recall an importat statement from matrix analysis 
LEMMA 1 For each pair of complex matrices A € C mx ™ anc / 
B € C" xm with m < n, the square matrix BA £ C" x ™ has 
the same eigenvalues as AB € ^yixm^ coun ti n g both alge- 
braic and geometrical multiplicity, together with an additional 
n — m eigenvalues equal to 0. 



g = r*r 



(4) 



G, 



(5) 



For the proof of this Lemma, we remand to |6|. Now, we 
prove the following 

THEOREM 1 The fco + ki non-zero eigenvalues rj k of the 
decision operator D = q\P\ — qapo can be evaluated as 
the eigenvalues of the SGM G s defined by (0. Moreover, the 
related eigenvectors |%) of D then the optimum measurement 
operators from (0, can be evaluated from the eigenvectors 9 k 
of the SGM for the same eigenvalue 7]k as \r/k) = c k T0 k where 
the factor c k — 1/ ^J6"^GQ k assures (rj k \rj k ) = 1. 

In words, the eigenvalues required to evaluate the correct 
detection probability according to ©, can be obtained by 
solving the eigenvalue equation 

—Goo — vlk —Goi 
Gio Gn — Jjifej 



det 







(6) 



where I k denotes the identity matrix of order k 

I k 



Proof: Let / 



By defining A = IT*, it is easy 
to prove that D = TA and G s = AT. Now, from the Lemma 
1, rjk, counted with their multiplicity, are eigenvalues of G s . 

To prove the second part note that \atk) = T9 k is an 
eigenvector of D. In fact D\ctk) = DT6 k = TG s 9k = 
VkFOk = r}k\otk)- Now we have to show that if |%) is a non- 
zero eigenvector of D it can be obtained by an appropriated 
6 k . Let 6 k = A\r) k ), we have that G s 6 k = G s A\ri k ) = 
AD\rj k ) = r] k A\ri k ) = r/ k 9 k then it is an eigenvector for 
G s and it cannot be zero, because otherwise A\r] k ) = => 
TAI^fc) = D\rj k ) — r] k \rj k ) = with contradiction. ■ 

Remark. We observe that if T is full-rank, then SGM is 
full-rank from rank(r) = rank(rT) = rank(G s ). In this 
case, rank(D) = rank(G s ) follows from Lemma 1. But it is 
not always true: in general we have lm(D) C Im(r), then 
rank(D) < rank(r) = rank(G s ). A simple example which 
leads to rank(£>) < rank(G s ) is provided by 
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with k = ki = 2. (7) 
3 whereas 



P'o 



Here T is not full-rank matrix and rank(G s ) 
rank(D) = 2. 

III. Check of the result with pure states 

In the case of rank 1 (pure states) the weighted density 
operators have the form 

q \a)(a\, p' x = qi \b){b\ (8) 

where \a) and \b) are normalized pure states. 

A. Helstrom's approach 

The decision operator D is given by 

D = p' 1 -p' Q = qi\b)(b\-q \a)(a\ (9) 

and has rank 2 (if \b)(b\ and \a)(a\ are linearly independent). 
For the EID of D we have to know the components of the 
kets | a) and \b), namely 

| a) = [a a\ ■■■ a„_i] T , \b) = [b bi ■■■ fo„_i] T . (10) 



Hence Dij — qib*bj — qoa*a,j, i,j = 0, 1, . . . , n — 1, where 
the components are linked by the normalization conditions 
(a\a) = (b\b) = 1 and qo+qi = 1. But the explicit evaluation 
of the eigenvalue equation det(D — rjl) is not immediate, 
since it involves a determinant of an n x n matrix. So, another 
approach is followed [4], based on the fact that D has rank 2 
and the eigenvectors 1?^) must belong to the two dimensional 
subspace spanned by D, so they are a linear combination 
of the states, namely \rji) = Ai\a) + Bi\b). Introducing this 
linear combination into the eigenvalue/eigenvector definition 
D\r]i) = rji\rji) and considering the normalization condition, 
one gets g 1 (A i X* + S i )|6) - q^A. +B, X)\a) = miMa) + 
Bi\b)) where X = (a\b). Finally, considering that the states 
are linearly independent, so that their coefficients must be 
equal, we obtain the equations 

q 1 {A i X*+B i ) =r H B i 
qo(Ai + Bi X) = rn Ai ' 

whose solutions give the non-zero eigenvalues: rj± = ^{qi 



i = 1,2 



(ID 



qo =F-\/l — 4qoqi\X \ 2 J, where rj + > 0, 77 < 0. Hence, from 
dij the correct decision probability is P c = qo + rj + , that is 



Pr = 



1 + Vi-^ogil^l 2 



which represents the so called Helstom 's bound. 
B. Skew Gram matrix (SGM) approach 

Considering the normalization, (a\a) = (b\b) 
notation X = (a\b), the skew Gram matrix is 

-qo -^qaqi X 
y/qoqi x* qi 

Then, equation (O becomes 



(12) 



1, and the 



G s 



det 



-qo ~ V 

y/q qiX* 



-■s/qooTX 
qi-v 



= 



(13) 



(14) 



which has exaclty the solutions t]± obtained above. 

Comparison. We have seen that Helstrom's approach is very 
articulated, because it starts with data expressed in terms of 
outer products, whereas the final result is expressed in terms 
of inner products (the parameter X = (a\b) in the case of 
pure states). The SGM approach is straightforward because it 
starts directly with inner products, the same we find in the final 
result. Note also that in the Helstrom's approach the original 
data are redundant (states of size n), whereas the system 
performance are completely determined by the geometry of 
inner products. 

IV. An application to quantum state comparison 

A problem of particular interest is quantum state compar- 
ison, where one wants to determine whether the states of a 
quantum system are identical or not 0. It can be formulated 
as a binary problem of discrimination between a pure state | a) 
and a uniformly mixed state, namely 

1 h 

p Q = \a){a\, Pl = -J2\h){bi\ (15) 

i=l 



where the h kets are supposed orthonormal, i. e. (bi\bj) = §ij. 

We evaluate the probability of correct discrimination using 
the SGM approach. The factors of the above density operators 
are 70 = y/qb~ \a) and 71 = y/qi/h [\h),. . . , \b h )} and the 
corresponding SGM is 



G s 



-qo -V 
V* (qi/h)I h 



(16) 



where V = ^Jq qi/h [X\, . . . , X h ] with X t = (a\h) and I h 
is the h x h identity matrix. The eigenvalue equation © gives 



qo) 



qoqi 



\x\ 



= (17) 



where ||X|| = J2i=i l^| 2 - ^ ts h + 1 solutions are 



qi ■ , 

m = -r , & = 1, • 



,h-l 



V± 



qi 
h 



qo ± 



qo 



Moqi 



\X\ 



(18) 



and are all positive, ry_ excepted. Hence, (0 gives 
h-1 



1 



1 + 91- 



91 , 



— 11*11 



(19) 



The case of interest is when all the h + 1 kets in ( fT5l ) have the 
same probability, that is when q\ — h/(h + 1). In this case 

l/(h- 



the error probability is P e = ll(h + 1) I 1 - Wl - l-Xl 

Comparison with J7J, where the direct Helstrom approach 
is used, shows that the SGM approach is much more simpler. 

V. General case of rank-2 density operators 

We develop the performance evaluation of the binary quan- 
tum system when the density operators po and p\ have rank 2. 
A particular attention is paid for the reduction of the number 
of parameters to get compact and readable results. 

A. Helstrom 's approach 

The weighted density operators have the form p' = 
qo {Pa \a)(a\ +Pc \c)(c\), p\ = qi (p b \b)(b\ +p d \d)(d\) where 
p a , pt, p c = 1 — p a , and pa = 1 — pb are the probabilities of 
the normalized states \a), \b), |c), and \d), respectively. In an 
rt-dimensional Hilbert space the states have n components: a;, 
bi, Ci, di, i = 0, 1, . . . , n — 1. Without loss of generality we 
can assume that the states forming the same density operator 
are orthogonal (8), that is 

(a\c)=0, (b\d)=0. (20) 

With the above specifications the decision operator becomes 

D = qipi - q p 

= qi [Pb \b)(b\+p d \d)(d\] - q [p a \a)(a\+p c |c)(c[ 

and its entries are = q x [p b bib* + p d did*} - q Q [p a aja* + 
p c CiC*]. Considering that qo + qi = 1, p a + Pc = 1 and 
Pb+Pd = 1, the specification of the decision operator is given 
by qo, p a , pb, I a), \b), |c), \d), where the kets are normalized 
and verify conditions (l20l . 



(21) 



For the evaluation of the eigenvalues of D we follow the 
procedure seen with pure states. Now, D has rank 4 and the 
eigenvectors |r?,) must belong to the four dimensional subspace 
spanned by D, so they are a linear combination of the states, 
namely \rji) = Ai\a) + Bi\b) + Ci\c) + Di\d). Hence, intro- 
ducing this linear combination into the eigenvalue/eigenvector 
definition D\rji) — r)i\r)i) and considering the normalization 
condition and orthogonality (l20l . one gets 

Ai qiPa \a) - Ai[q p b X* \b) + Pd Y* \d)]+ 
+B iqi [p a X\a) + Pc W\c)}- B t qop b \b) + 

+a qi Pc \ c ) - a qo[ Pb w* \b) + Pd z* \d)}+ (22) 

+Di qi[p a Y I a) + p c Z \c)] - Di q p d \d) 
=r)i[Ai I a) + B t \b) + C l |c) + D t \d) 
where we have introduced the inner products 
X = (a\b), Y = (a\d), W = (c\b), Z = (c\d) . (23) 

Finally, considering that the 4 states are linearly indepen- 
dent, so that their coefficients must be equal, we obtain 4 
equations in the unknowns Ai, Bi,Ci, Di and rji. By solving 
with respect to % we obtain an algebraic equation of degree 4 
whose solutions give the non zero eigenvalues. The procedure 
is cumbersome and we do not write explicitly the passages, 
since they are immediate from the SGM procedure. 

B. Skew Gram matrix approach 

We have to find the factors of the density operators. Con- 
sidering the orthogonality condition (l20l . they are given by 

7o = qo [VK\a), VPc~\c)} , 71 = 1i [y/Pb \b), y/pd\d)] ■ (24) 

Considering again orthogonality ( 1201 ) and normalization, the 
blocks of the SGM are given by 

Pa 
Pc 



Qo 



Goo = 7o7o 

Gn = 7*7i = 1i 
where 



x 
w 



Pb 






Pd 



G01 = 7o7i 
Gio=7i*7o = Qi, (25) 



x 
w 



\/qoqiPaPbX, 



yJqoqiPcPbW, 



^/qoqiPaPdY, 



y/qoQlPcPdZ, 



(26) 



are weighted inner products. We see that in general the SGM 
G s depends on the seven parameters qo, Pa , pb, x, y, w, z. 
The eigenvalue equation of G s is 

~ 4 1 Brf + Cr) 2 + Drj + E = 



V 



(27) 



where 



\y\ 2 ps 



E = \xz\ + \wy\ — \x\ qs — \w\ qr - 

— xw* zy* — x*wz*y + pqrs 
D = \x\ 2 (q -s) + \w\ 2 (q - r) + \y\ 2 (p - s)+ 

+ \z\ 2 (p-r)- 
C=\x\ 2 + \w\ 2 + \y 
B=p + q-r-s = q -qi 



\z\ 2 pr+ 



- prq + prs 

2 + \z\ 2 +pq—pr — qr—ps — qs + rs 



(28) 



with 

p = qoPa, q = qoPc, r = qm, s = q x p c . (29) 

Equation j27| ) is a quartic equation, whose analytical so- 
lutions are known from the times of Cardano, but their 
expression are very long. When the states have equal prior 
probabilities we find B = and d27T i becomes the depressed 
quartic equation for the absence the third degree term. The 
depressed quartic equation played an important role in the 
history of mathematics. However, we need the sum 7/1 + 7/2 
of the two positive solutions, which is more readable than the 
individual solutions. We find 



m + 7 n = 



\ 



~~~ + 3?2 



^2S 



(30) 



B 



3t/R+VR 2 -4S 3 2 
27D 2 + 21B 2 E and 



where R = 2G 3 - 9BDC - 72EC 
S = C 2 - WD + 12E. 

C. Interpretation of parameters and numerical example 

We can take as reference the case of pure states \a) and 
I b) and their inner product X = (a\b), which determines the 
error probability through \X\ 2 according to (112t . In the ideal 
case the states are orthogonal, X = 0, and give P c = 1. When 
\X\ > we find P c < 1 and the degradation is to ascribe to 
the presence of the shot noise (in the classical interpretation). 

Now, it is difficult to choose the data listed in d29l to 
carry out a numerical example, owing to the several conditions 
therein, which ultimately assured that the correspondent Gram 
matrix is PSD. To be sure that the data are feasible, we 
have chosen the density operators from Glauber's theory on 
coherent states, modeling the monochromatic electromagnetic 
radiations produced by a laser. In these theory a density 
operator p(a), represented by a matrix of infinite dimension, 
depends only on two parameters a and J\f, with N a = \a\ 2 
giving the average number of signal photons, and JV the 
average number of thermal noise photons. For Af = (absence 
of thermal noise) we find rank 1, that is pure states. Thus, to 
find rank 2 we have to consider a small amount of thermal 
noise. Note that for practical calculations, we need a finite 
n-dimensional approximation, of the infinite dimensional rep- 
resentation. As discussed in detail in [9|, to get a moderately 
small value of the size n we have to choose a small value of 
N a = \ct\ 2 and a rank 2 is assured with a very small value 
of J\f. In practice, a good choice for a non symmetric case 
may be n = 10, q = 0.4, a Q = -1.2247 -> N ao = 1.5, 
«! = 1.3038 -> N ai = 1.7, TV = 0.05. 

With these parameters we have evaluated the 10 x 10 density 
operators po = p(ao), p\ — p(a±), and the 10 x 2 factors 70, 
71 . For reason of space we omit the explicit numerical matrices 
and we write directly the blocks of the SGM, which result 



Go 



G11 — 



0.381 


0.571 






0.018 


0.027 



G01 = 



Gin — G, 



0.019 
0.011 

01 • 



0.011 
0.005 



From OH) we obtain (see J25l) ) the probabilities p a = 
0.95240, p c = 0.04534, p b = 0.95243, p d = 0.04533 and the 
inner products X = 0.04089, Y = 0.10346, W = 0.10345, 
Z = 0.22043. We have evaluated the 10 eigenvalues of 
D, which result {-0.380307, -0.0174054, 0, 0, 0, 0, 0, 0, 
0.0263829, 0.570885} and give P c = 0.997268 , P e = 
1 - P c = 0.00273197. We have checked that the 4 eigen- 
values of G s are exactly the non zero eigenvalues of D, in 
agreement with the theory. The coefficients of the eigenvalue 
equation ([27]> are E = 0.000105191, D = 0.00183749, 
G = -0.216064, B = -0.198617 and the solutions are 
{-0.380307, -0.0183563, 0.0263943, 0.570885}. To get the 
sum of the positive solutions d30T > from the formulas, we 
have evaluated R = -0.0190434, S = 0.0490408. Then 
V1+V2 = 0.59728, which gives again tha above probabilities. 

The error probability with pure states, obtained with the 
same inner product X = 0.04089, is P e = 0.000401349, 
that is one order of magnitude better than with rank 2. The 
degradation with mixed states is due to the presence of thermal 
noise. 

VI. Rank 2 density operators with symmetry 

An m-ary state constellation {70, 71, • • • ,7m-i} exhibits 
the geometrically uniform symmetry (GUS), see Q5J , if the 
states are related as 7$ = S* l 7o, where S is a unitary operator 
S, such that S m = I-u is the identity operator in the Hilbert 
space %. With the GUS equal prior probability are assumed, 
that is qi = 1/m, 

Helstrom's approach is not very much simplified by the 
GUS, so we pass directly to the alternative equivalent ap- 
proach. 

A. SGM approach with GUS 

We investigate the simplifications on the Gram matrices due 
to the GUS and equal prior probabilities. We know that, in the 
presence of GUS, the Gram matrix becomes block-circulant, 
that is its blocks Gy depends only on the index differences 
j — i |9|. In the binary case the block-circulant conditions are 
G11 = Goo, G01 = Gio = Gq X . Inspection on (l25l l shows that 
this condition leads to the following simplifications: I) Pb = 
p a and pd = Pc, 2) the inner product X = (a\b) , Z = (c\d) 
are real, 3) the inner products Y and W are conjugate. Hence, 
the leading blocks take the form 



Goo — 7: 



Pa 




G01 — - 



PaX 
^/PaFcY* 



VPaPc Y 
PcZ 



(32) 





Pc 

where p a + p c = 1. In conclusion with the GUS the SGM 
depends only on the four parameters p ai X, Y, Z. 

The eigenvalue equation of G s , given by ( fl4] l. written in a 
convenient form, is 



±L = Q 



(31) 



V A -\H V 2 
where 

H = p 2 a (l - X 2 )+p 2 c (l - Z 2 ) - 2p aP c\Y\ 2 
L = ( PaPc ) 2 [\Y\ 4 - 2(1 + XZ)\Y\ 2 + (1 - X 2 ){1 
= (p aPc ) 2 [(\Y\ 2 - (1 + XZ)) 2 -(X + Z) 2 ] . 



(33) 



Z 2 )] 
(34) 



Hence, the quartic equation d27b degenerates into a bi- 
quadratic equation, whose solutions are straightforward. In 
particular, the two positive solutions are given by T]± t 2 = 

\{H± VH 2 - AL) , hence 



Pn 



2 ' 2 V 2 



H + 2VZ. 
The explicit result is 



UH + s/H 2 -4L) + y±(H - VH 2 - 4L) 



(35) 



Pc = \ + \ (p 2 a (l-X 2 )+pl(l-Z 2 )-2p aPc \Y\ 

+2 P aPc y/[\Y\ 2 -(l+XZ)]*-{X + Z)2 



1/2 



(36) 



Check with pure states. The case of pure states is obtained 
by letting p a = 1 and then p c = 0. In agreement with {12}, 
from (l36l l we get 



Pr 



(37) 



Case of orthogonality. If Y = (a\d) = 0, also W = (c\b) = 
0. This leads to a simplification in d36l l. namely 



P = i 



aVl-X 2 + PcA /l- Z 2 



(38) 



B. Numerical example 

A choice for a symmetric case may be n = 10, qa = 0.5, 
a Q = -1.26491 -> N ao = 1.6, ai = 1.26491 -> 7V ai = 1.6, 
TV = 0.05 then, the leading blocks of the Gram matrices are 



Goo — 



7o7o 



Got = 7q7i 



0.476206 
0.0226692 

0.019409 -0.0107206 
-0.0107206 0.00498695 



(39) 



From 
0.95241 
and the 

W = 



we obtain (see (l25l l) the probabilities: p a — 
p c = 0.04534, p b = 0.95241, p d = 0.04534 
inner products X = 0.04076, Y = -0.10318, 
-0.10318, Z = 0.21999. The eigenvalues of D 
are {-0.47558, -0.0218743, 0, 0, 0, 0, 0, 0, 0.0218743, 
0.47558} and give P c = 0.997268 and P e = 0.00273197. 
The eigenvalues of G s are exactly the non zero eigenvalues 
of D, in agreement with the theory. The coefficients of 
the eigenvalue equation (03) are -j$L = 0.000108222, 
\H = -0.226655 and the solutions {-0.47558, -0.0218743, 
0.0218743, 0.47558}. 

We now check that the same results are obtained from 
the analytical formulas. The evaluation of (l34b gives H = 
0.907538 and L = 0.0017714. Hence, from ([36) we obtain 
P c = 0.997924, in agreement with the above evaluation. 

The case Y = gives H = 0.907538, L = 0.0017714, and 
P c = 0.997924, P e = 0.00207581. 

The error probability with pure states (with the same inner 
product X = 0.04076) is P e = 0.000415467, that is one 



order of magnitude better than with with with rank 2. The 
degradation is due to the presence of thermal noise. 

C. Limit to closed-form results 

As established by Evarist Galois two hundred years ago, the 
solution of algebraic equations can be written explicitly up to 
the fourth order. We have seen that in the general case of the 
binary quantum detection the case of rank 2+2 leads just to a 
quartic equation and so it is the limit. In the presence of GUS 
the quartic equation is essentially reduced to the second order. 
In principle, with the GUS it would be possible to solve the 
case of rank 3+3, where the six order equation can be reduced 
to a cubic, and also the case of rank 4+4, where the eight order 
equation can be reduced to a quartic. No more closed form is 
possible, unless in particular cases. 

VII. Conclusions 

The Helstrom theory on binary detection is the only general 
explicit result available in quantum optimization, but its trans- 
lation into formulas is really possible only with very small 
dimensions of the density operators. On the other hand, the 
final results depends only on the inner products of quantum 
states, which compress the information contained in the den- 
sity operators. The skew Gram matrix approach, which gives 
exactly the same results, starts just from the inner products, 
with a dramatic simplification of the algebra involved, and 
permits to establish closed-form results not available else- 
where. Furthermore, the skew Gram matrix approach clearly 
establishes a fundamental truth, at least for binary detection: 
the inner products, collected in the Gram matrix, give the 
necessary and sufficient information to achieve the optimal 
detection. 
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